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Abstract 

Challenges  associated  with  the  allocation  of  limited  resources  to  mitigate  the  impact  of  natural  disasters  inspire 
fundamentally  new  theoretical  questions  for  dynamic  decision  making  in  coupled  human  and  natural  systems.  Wildfires  are 
one  of  several  types  of  disaster  phenomena,  including  oil  spills  and  disease  epidemics,  where  (1)  the  disaster  evolves  on  the 
same  timescale  as  the  response  effort,  and  (2)  delays  in  response  can  lead  to  increased  disaster  severity  and  thus  greater 
demand  for  resources.  We  introduce  a  minimal  stochastic  process  to  represent  wildfire  progression  that  nonetheless 
accurately  captures  the  heavy  tailed  statistical  distribution  of  fire  sizes  observed  in  nature.  We  then  couple  this  model  for  fire 
spread  to  a  series  of  response  models  that  isolate  fundamental  tradeoffs  both  in  the  strength  and  timing  of  response  and 
also  in  division  of  limited  resources  across  multiple  competing  suppression  efforts.  Using  this  framework,  we  compute 
optimal  strategies  for  decision  making  scenarios  that  arise  in  fire  response  policy. 
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Introduction 

Wildfire  progression  is  driven  by  inherently  complex,  inter¬ 
connected,  physical  processes,  involving  a  variety  of  factors, 
including  weather,  vegetation,  and  terrain  [1—4],  Wildfires  are 
difficult  to  predict  in  space  and  time,  and  thus  it  is  hard  to  make 
informed  decisions  about  mitigation  efforts,  particularly  the 
allocation  of  limited  resources  during  dynamic  response.  Further¬ 
more,  the  areas  and  damages  of  wildfires  both  exhibit  high 
variability,  heavy  tailed  statistics  [5] ,  implying  that  overall  impacts 
are  dominated  not  by  the  typical,  median  size  events,  but  rather  by 
the  relatively  rare,  large  events.  This  highlights  the  fundamental 
importance  of  understanding  tradeoffs  between  strength  and 
timing  in  mitigation  efforts  that  may  help  prevent  future  wildfires 
from  growing  out  of  control. 

In  recent  years,  wildfires  in  the  United  States  have  received 
media  attention  due  to  rising  suppression  costs  [6],  which  have 
regularly  exceeded  expenditures  of  a  billion  dollars  per  year  in  the 
new  millennium  [7],  A  variety  of  potential  driving  factors  have 
been  suggested  [8],  including  a  climate-driven  increase  in  wildfire 
activity  [9-11],  and  development  at  the  wildland  urban  interface 
[12,13],  A  cohesive  strategy  that  takes  into  account  both  the 
natural  occurrence  of  wildfire  on  the  landscape  and  its  impact  on 
human  lives  is  needed  to  determine  the  best  course  of  action  for 
reducing  these  costs.  The  Federal  Land  Assistance,  Management 
and  Enhancement  Act  (FLAME  Act)  of  2010  directs  the  U.S. 
Department  of  Agriculture,  the  U.S.  Department  of  the  Interior, 
and  the  U.S.  Department  of  Homeland  Security  to  develop 
a  Cohesive  Wildfire  Management  Strategy  to  identify,  among 


other  things,  the  most  cost  effective  means  for  allocating  fire 
management  budget  resources.  Thus,  there  is  a  need  to  un¬ 
derstand  not  only  the  dynamics  underlying  wildfire  behavior  but 
also  how  to  allocate  limited  resources  to  mitigate  fire  damage. 

In  this  paper  we  introduce  a  minimal  model  of  wildfire 
dynamics  using  the  formalism  of  stochastic  processes  and  queueing 
theory.  We  divide  the  burnable  substrate,  which  may  be  forest, 
shrubland,  or  grassland  depending  on  the  region,  into  discrete 
parcels  that  arrive  in  a  queue  when  they  begin  to  burn.  The 
growth  rate  of  the  number  of  parcels  in  the  queue  scales  with  the 
current  size  of  the  queue,  capturing  the  escalating  dynamics 
familiar  in  wildfire  spread.  Parcels  exit  the  queue  when  they  are 
extinguished  either  by  natural  causes  or  human  suppression.  This 
model  reproduces  the  observed  statistical  distribution  of  wildfire 
sizes,  and  it  enables  investigation  of  tradeoffs  in  suppression  efforts 
that  are  fundamental  to  decision  making. 

Our  model  is  sufficiently  general  that  it  may  apply  to  other 
scenarios  involving  coupled  dynamics  of  a  disaster  and  response. 
In  the  context  of  wildfires,  our  study  is  intended  to  be  an 
intermediate  step  between  models  with  static  fires  [14,15]  and 
models  with  a  detailed  fire  spread  simulation  [16,17],  It  is  not 
intended  to  replace  high  fidelity  simulations  but  will  be  used  to 
examine  basic  relationships,  which  can  then  be  further  explored 
and  verified  with  more  realistic  models.  We  comment  on  this  in 
more  detail  in  the  Discussion. 

We  focus  specifically  on  time  dynamics.  However,  spatial  factors 
such  as  prepositioning  of  resources  are  still  implicitly  included 
through  the  time  delay  of  resource  arrival.  Tradeoffs  in  strength 
and  timing  are  particularly  relevant  when  considering  resource 
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allocations  to  simultaneously  occurring  fires,  which  may  be 
evolving  on  different  time  scales. 

Our  Contributions  in  Context 

There  are  many  research  efforts  that  are  complementary  to 
ours,  having  focused  primarily  on  spatial  rather  than  temporal 
aspects  of  the  problem.  A  rich  literature  exists  within  operations 
research,  focusing  on  optimization  of  wildfire  response  [18], 
including  models  with  realistic  fire  spread  and  various  types  of 
suppression  resources  [16,17,19],  Because  of  the  computational 
complexity  inherent  in  simulating  a  single  fire,  much  less  multiple 
co-occurring  fires  or  a  statistical  distribution  of  fires  over  time, 
these  models  often  abstract  away  the  time  dynamics  of  individual 
fires.  In  many  existing  studies,  the  decision  variable  relates  to 
spatial  prepositioning  of  resources  for  the  initial  response. 
However,  because  wildfires  are  fundamentally  both  geospatial 
and  dynamic,  decision  making  tools  must  ultimately  capture  both 
the  spatio-temporal  spread  of  the  fire  itself  and  the  geographic 
locations  and  transportation  of  suppression  resources. 

One  aspect  of  wildfires  that  has  received  considerable  attention 
in  the  complex  systems  literature  is  the  power  law  distribution  of 
fire  sizes:  P(F)  where  F is  the  footprint  (i.e.,  total  burn  area) 

of  the  fire  and  P(F)  is  the  probability  of  a  fire  being  larger  than  or 
equal  to  size  F.  Fires  are  not  unique  in  this  regard;  many  disparate 
disaster  phenomena  exhibit  power  law  statistics  when  losses  are 
tallied  [20,21],  and  many  biological,  technological,  and  social 
systems  also  exhibit  power  laws  [22-24],  Although  these 
observations  have  led  to  widespread  association  of  power  law 
statistics  with  a  key  signature  of  self-organization  in  complex 
systems  [25,26],  there  are  mathematical,  statistical,  and  data- 
analytic  arguments  that  suggest  power  law  distributions  should  be 
the  natural  null  hypothesis  for  high  variability  phenomena  [27]. 

Power  law  statistics  are  by  no  means  unique  or  unexpected  for 
wildfires.  However,  the  impact  of  power  law  distributions  for 
policy  and  decision  making  has  not  been  investigated  in  much 
detail,  yet  is  critically  important,  especially  in  the  case  of  natural 
disasters.  This  heavy  tailed  statistical  form  is  intrinsically  linked 
with  uncertainty  and  high  variability  in  observations:  most  fires  are 
small,  while  most  of  the  area  burned  is  contained  in  the  few  largest 
events.  In  fact,  in  the  context  of  wildfires,  it  has  been  noted  that  the 
top  1%  of  large  fires  burn  80—96%  of  the  total  area  [28],  Another 
consequence  of  this  mathematical  form  is  that  a  large  event,  many 
times  bigger  than  anything  seen  recently,  is  not  only  possible  but  is 
also  consistent  with  the  distribution.  This  is  particularly  significant 
when  planning  for  a  worst  case  scenario. 

Although  power  law  statistics  are  observed  in  many  systems, 
a  particularly  interesting  feature  is  that  the  exponent  for  wildfires, 
a,  clusters  around  1/2  for  data  sets  spanning  a  wide  range  of 
geographic  regions  and  time  frames  [4],  In  general,  different 
disaster  types  (e.g.,  fires,  hurricanes,  earthquakes)  exhibit  varying 
values  of  the  exponent,  suggesting  that  the  relative  scaling  of 
physical  processes,  tradeoffs,  and  constraints  may  be  different  in 
each  case  [29—31].  The  consistent  value  of  the  exponent  a  for 
wildfire  distributions,  despite  varying  vegetation,  weather,  and 
terrain,  suggests  that  there  may  be  a  common  mechanism  that 
drives  wildfires  across  a  wide  range  of  fire  regimes. 

Several  mechanisms  for  power  laws  have  been  discussed  in  the 
context  of  wildfires,  and  simple  models  have  been  introduced  to 
illustrate  these  mechanisms.  One  example  is  Self  Organized 
Criticality  (SOC)  [32,33],  which  aims  to  explain  the  presence  of 
power  laws  as  the  consequence  of  a  dynamical  system  operating  at 
a  critical  phase  transition.  Sandpile  cellular  automata  are  the 
canonical  example  of  SOC  [25],  combining  an  infinitesimal  rate 
for  discrete  and  spatially  random  reloading  of  the  system  with 


a  threshold  rule  for  event  propagation,  resulting  in  power  law 
distributions.  Extrapolation  of  this  basic  mechanism  of  slow, 
random  regrowth  with  fast  event  propagation,  terminating  at  grid 
locations  where  the  threshold  criterion  for  propagation  is  not  met, 
has  made  SOC  a  popular  topic  within  the  complex  systems 
community. 

The  phase  transition  underlying  the  SOC  forest  fire  model  is 
a  percolation  transition,  associated  with  a  critical  density  at  which 
nearest  neighbor  connectivity  first  extends  from  one  side  of 
a  discrete  lattice  to  another.  The  percolation  transition  is 
fundamentally  statistical,  describing  the  average  behavior  of  an 
ensemble  consisting  of  all  possible  configurations  of  occupied  (by 
fuel)  and  vacant  (no  fuel)  sites,  characterized  only  by  the  fuel 
density.  Percolation  is  the  most  widely  studied  model  of  criticality 
in  statistical  and  mathematical  physics  because  it  is  so  simple, 
enabling  rigorous  analysis  and  large  scale  simulations  [34] .  In  the 
context  of  wildfires,  the  SOC  model  can  be  loosely  interpreted  as 
a  fuel  driven  model,  because  fire  propagates  through  connected 
clusters  of  fuel  occupied  sites,  terminating  only  at  vacancies. 
However,  even  in  a  physical  model  of  fuel  driven  wildfires,  the 
biology,  ecology,  chemistry,  and  physics  of  fuels  is  more  complex 
than  what  is  captured  by  an  ensemble  of  random  configurations  of 
occupied  and  vacant  sites.  It  is  therefore  no  surprise  that  the 
simple,  mechanistically  illustrative  SOC  forest  fire  model  fails  to 
describe  wildfire  statistics  (the  exponent  does  not  agree  with 
observations)  [35],  Taken  literally,  the  SOC  model  would  also 
predict  wildland  vegetation  densities  corresponding  to  the  critical 
density  in  percolation,  which  is  much  lower  than  realistic 
measurements  for  burnable  terrain  [36,37],  as  well  as  fire 
footprints  that  are  tenuous  and  fractal,  as  opposed  to  the  observed 
compact  fire  shapes  [38] .  Thus,  the  system  would  remain  poised  at 
the  critical  density  even  immediately  after  a  large  region  spanning 
fire. 

More  recently,  Highly  Optimized  Tolerance  (HOT)  [35,39] 
was  introduced  to  illustrate  characteristics  generally  expected  to 
arise  in  optimized  (e.g.,  technology),  organized  (e.g.,  ecology),  or 
evolved  (e.g.,  biology)  systems  subject  to  resource  constraints  and 
a  spectrum  of  external  disturbances.  The  surprising  result  is  that 
this  abstract  problem  definition  alone  also  naturally  leads  to 
a  mechanism  for  power  laws.  The  HOT  mechanism  was  described 
initially  using  fire  spread  terminology,  in  order  to  clarify  the 
differences  between  HOT  and  SOC.  The  simplest,  most  abstract 
HOT  wildfire  model  posits  that  power  laws  emerge  due  to 
tradeoffs  between  vegetation  yield,  cost  of  resources,  and  tolerance 
to  risks,  whether  that  tolerance  is  a  compromise  made  by 
corporations  in  a  managed  forest,  or  whether  it  involves  biological 
and  ecological  tradeoffs  between  growth  and  resilience,  adapted  to 
conditions  in  a  particular  fire  prone  terrestrial  ecosystem.  Losses 
are  represented  as  “area  burned,”  and  “resources”  represent  the 
effort  to  minimize  average  losses.  This  model  thus  does  not  assume 
that  fuel  availability  is  the  only  determining  factor  for  fire  size 
distributions,  but  rather  that  it  emerges  as  a  consequence  of  some 
organizational  process  that  over  time  favors  maximization  of  the 
average  yield  or  proliferation  of  vegetation  well-suited  to  the 
regime.  Compared  to  SOC,  HOT  is  much  more  efficient  in 
utilization  of  resources.  In  the  context  of  wildfires,  simple  HOT 
models  produce  compact  two-dimensional  fires,  extinguishing  at 
the  one-dimensional  perimeters.  HOT  leads  to  high  vegetation 
densities,  much  more  representative  of  natural  wild  lands,  and 
HOT  also  leads  to  heavy  tailed  size  distributions.  This  can  be 
understood  because  optimizing  the  system  for  common  events 
leads  to  catastrophically  large  rare  events.  Interestingly,  HOT 
produces  a  power  law  that  matches  both  wildfire  data  (a  =1/2) 
and  detailed  fire  simulations  [4,40-42]. 
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In  this  work  we  posit  an  alternate  simple  model  in  which  we 
focus  on  temporal  dynamics  rather  than  spatial  features.  Our  first 
objective  is  to  understand  the  extent  to  which  the  dynamics  of  fire 
growth  and  suppression  can  be  understood  from  models  that 
capture  the  temporal  features  but  abstract  the  spatial  features.  The 
second  objective  is  to  use  this  simplified  temporal  framework  to 
assess  tradeoffs  that  result  from  basic  tensions  in  fire  response 
policy:  fire  growth  over  time,  limited  resources,  and  decisions 
about  when  to  deploy  resources  and  how  many  to  send.  Our  study 
takes  an  important  preliminary  step  in  the  development  of  future 
applications  that  facilitate  real-time  response  efforts  in  coupled 
human  and  natural  systems. 

Methods 

Natural  Fire  Dynamics 

We  begin  with  a  minimal  model  of  fire  growth  and  natural  fire 
extinction,  based  on  the  physical  processes  known  to  drive  fire 
behavior.  In  the  next  section,  we  augment  the  model  to  include 
human  suppression  explicitly.  To  capture  the  temporal  dynamics 
of  wildfire,  we  consider  a  finite  number  of  discrete  parcels  of 
burnable  substrate  that  can  be  affected  by  fire.  We  associate  each 
parcel  with  one  of  three  states:  unhurried,  burning,  and  burned.  Our 
simulations  start  with  some  number  of  burning  parcels  and  the  rest 
of  the  parcels  in  the  unburned  state.  An  unburned  parcel  can  catch 
fire  to  enter  the  burning  state  and  will  eventually  die  out, 
becoming  burned.  The  time  scale  under  consideration  here  does 
not  allow  for  regrowth  of  a  burned  parcel,  so  these  state  transitions 
proceed  in  one  direction  only. 

We  refer  to  burning  parcels  as  firelets  and  assume  that  they  can 
spread  to  unburned  parcels,  i.e.  firelets  can  give  “birth”  to 
additional  firelets.  Although  the  notion  of  discrete  parcels  is 
geospatially  motivated,  our  rules  for  fire  spread  ignore  spatial 
features.  However,  spatial  scale  is  implicit  in  our  model  in  two 
important  ways.  First,  the  rate  at  which  a  firelet  gives  birth  to  other 
firelets  assumes  an  underlying  length  scale  associated  with  the  size 
of  the  parcel  and  the  rate  at  which  fire  spreads  across  it.  Second, 
the  total  number  of  parcels  J,  which  we  refer  to  loosely  as  the  size  of 
the  burnable  region,  sets  an  upper  bound  F<J  on  fire  footprint  F 
represented  in  our  model.  Note  that  J  need  not  be  interpreted 
literally  as  the  total  area  of  a  particular  region,  but  instead  may  be 
more  directly  linked  to  some  user  specified  threshold  area,  perhaps 
associated  with  a  scale  at  which  the  fire  exceeds  the  capacity  of 
locally  available  suppression  resources,  and  is  deemed  out  of  control. 
For  this  interpretation,  a  real  wildfire  may  continue  to  expand, 
consuming  an  even  larger  area,  but  growth  beyond  this  size  is  not 
explicidy  represented  numerically. 

A  finite  population  of  parcels  J  constrains  the  behavior  of  the 
model  in  two  ways.  First,  we  have  an  upper  bound  on  the  number 
of  potential  firelets.  As  the  number  of  unburned  parcels  decreases 
over  time,  the  upper  bound  on  firelets  also  decreases.  Second,  all 
fires  are  finite  in  duration  and  ultimately  end  with  the  fire  going 
out  (a  trapping  state),  because  the  fire  has  consumed  the  entire 
burnable  region  ( F  =  J )  or  because  the  fire  stopped  short  of 
burning  everything  ( F <J ).  In  either  case,  system  behavior  is 
transient,  so  steady  state  calculations  of  fire  behavior  are  not 
meaningful.  Instead,  we  examine  the  effects  of  model  parameters 
on  the  number  of  burned  parcels  F. 

We  model  the  transitions  for  a  population  of  parcels  as  state- 
dependent  and  stochastic.  Let  j  denote  the  current  number  of 
firelets.  Let  [l  be  the  constant  rate  at  which  a  single  firelet  gives  rise 
(or  birth)  to  a  new  firelet.  Let  (5  be  the  constant  rate  at  which  a  single 
firelet  is  extinguished  (or  dies)  on  its  own  (i.e.,  in  the  absence  of 
human  suppression).  Then,  the  state  transitions  of  the  population 


as  a  whole  can  be  understood  in  terms  of  a  birth-death  process 
where  in  state  j  the  aggregate  birth  rate  is  Ay  =  jf 1  and  the  aggregate 
death  rate  is  m  =  j5.  We  choose  the  simplest  assumptions:  all 
burning  firelets  are  actively  spreading  and  can  also  die  at  any  time. 
To  simulate  continuous  time  dynamics  we  model  the  time  interval 
until  the  next  birth  or  death  as  an  exponentially  distributed 
random  variable,  t.  The  probability  density  function  for  t  is: 
/(f)  =  (Ay  +  /t,)e_^+/‘y)(  This  implies  that  the  existence  of  more 
firelets  increases  the  probability  that  in  any  fixed  time  interval  at 
least  one  of  them  either  spreads  or  dies  out.  The  specific  choice  of 
rates  underlies  quantitative  relationships  (on  average)  between 
duration  and  size  F  and  the  characteristic  time  scale  associated 
with  system-wide  events. 

Time  Series  of  Individual  Fires 

Fire  behavior  is  driven  primarily  by  the  ratio  of  the  state- 
dependent  firelet  arrival  (Ay)  and  departure  (jl)  rates.  When 
Ay  <  Hj  the  fire  tends  to  decrease  in  size,  and  when  Ay  >  ft )  the  fire 
tends  to  increase  in  size.  Here  by  construction,  the  ratio  Ay/m  is 
constant  for  all  values  ofj  and  characterized  by  the  scaled  growth 
rate  [1/5,  which  serves  as  a  proxy  for  physical  conditions  such  as 
weather,  vegetation  moisture  content,  and  atmospheric  humidity. 
For  example,  high  winds  may  be  represented  by  an  increase  in 
the  birth  rate  /I,  and  a  higher  scaled  growth  rate.  Similarly,  high 
vegetation  moisture  may  be  represented  by  an  increase  in  the 
natural  death  rate  5,  and  a  lower  scaled  growth  rate.  For 
simplicity  we  drop  the  term  ‘scaled’  and  refer  to  [S/d  as  the  growth 
rate. 

In  the  analysis  to  follow,  we  select  growth  rate  values  that  yield 
size  distributions  roughly  matching  size  distribution  data.  In 
principle,  the  growth  rate  itself  could  be  estimated  for  a  given  set  of 
environmental  factors,  by  combining  empirical  measurements 
with  a  theoretical  understanding  of  fire  propagation.  For  example, 
the  Rothermel  equation  [43,44]  takes  wind,  vegetation,  and 
terrain  as  inputs  to  estimate  the  fire  spread  rate,  which  in  our  case 
corresponds  to  f).  Fuel  composition,  and  therefore  the  energy 
released  via  combustion,  would  determine  the  rate  <5  at  which  the 
fire  ceases  to  be  actively  spreading. 

Figure  1  illustrates  trajectories  for  two  simulated  fires.  At  each 
discrete  time  step  we  display  the  number  of  parcels  that  are  in  the 
burning  (i.e.,  they  firelets),  burned,  and  unburned  states.  As  one 
would  expect,  j  fluctuates  over  time,  with  the  fire  ending  when 
j  =  0.  Because  the  total  population  J  is  constant  and  finite,  the 
number  of  burned  parcels  is  nondecreasing  in  time  while  the 
number  of  unburned  parcels  is  nonincreasing. 

Figure  1A  shows  a  simulated  fire  with  growth  rate  [S/5<  1, 
a  condition  we  refer  to  as  mild.  This  fire  dies  out  with  a  nonzero 
number  of  unburned  parcels  remaining  ( F  <  J),  which  is  typical 
for  mild  conditions.  Figure  IB  shows  a  fire  with  growth  rate 
P/5>\,  which  we  refer  to  as  extreme  conditions.  The  simulated 
fire  continues  until  there  are  no  remaining  unburned  parcels 
(F  =  J).  Once  this  size  limit  is  reached,  the  number  of  firelets 
steadily  decreases  until  the  fire  is  out.  Fires  that  grow  out  of  control 
and  span  the  system  are  common  in  extreme  conditions.  In  both 
cases,  the  number  of  firelets  at  any  point  in  time  is  only  a  small 
fraction  of  the  total  number  of  burned  parcels,  of  order  j  ~  s/F, 
consistent  with  the  perimeter  of  a  compact,  two-dimensional  fire 
footprint. 

Results 

Fire  Size  Distributions 

Size  statistics  P(F),  obtained  by  simulating  catalogs  of  104 
individual  fires,  are  illustrated  for  different  growth  rates  in 


PLoS  ONE  |  www.plosone.org 


3 


April  2012  |  Volume  7  |  Issue  4  |  e33285 


A  Dynamic  Model  for  Wildfire  Suppression 


1000 

800 

<1> 

600 
z  400 
200 
0 


Figure  1.  Individual  fire  simulations  in  the  absence  of 
suppression  for  (A)  mild,  and  (B)  extreme  conditions.  The  total 
number  of  parcels  .7=1000  is  divided  between  unburned  (green), 
burning  (red),  and  burned  (black)  according  to  the  stochastic  dynamics 
of  the  birth-death  process.  The  fire  terminates  when  the  number  of 
burning  firelets  j  first  reaches  zero.  In  (A)  the  final  fire  size  F  is  roughly 
27/5,  while  in  (B)  F  =  J. 
doi:10.1371/journal.pone.0033285.g001 


Figure  2A.  For  mild  conditions,  fires  generally  go  out  before 
consuming  the  entire  burnable  region  (total  number  of  parcels 
7=1000),  and  P(F)  can  be  approximated  by  a  power  law  with 
a  large  scale  cutoff  characterized  by  some  Fmax<J.  In  contrast, 
for  extreme  conditions  fires  are  generally  much  larger  and  there  is 
an  accumulation  of  top  ranked  fires  that  burn  all  7  parcels.  This 
value  can  be  understood  as  a  threshold  beyond  which  the  fires 
cannot  be  controlled. 

Interestingly,  for  /S/<5  =  1  (a  case  we  denote  as  dangerous 
conditions),  P(F)  is  an  exact  power  law  with  exponent  1/2, 
similar  to  observed  wildfire  statistics  [4],  Mathematically,  the 
origin  of  this  exponent  can  be  understood  in  terms  of  a  discrete 
random  walk  in  which  each  “step”  corresponds  to  a  state 
transition  in  the  birth-death  process.  Our  P(F)  maps  directly 
onto  the  distribution  of  first  return  times  (the  number  of  firelets 
corresponds  to  displacement  of  the  random  walker,  and  fire  size 
maps  to  time),  which  is  a  power  law  with  exponent  1/2.  For 
a  mathematical  proof,  see  [45]. 

To  what  extent  does  our  model  capture  some  basic  features  of 
fire  spread  even  without  explicit  spatial  rules?  If  we  think  of  the 
number  of  firelets  as  the  region  of  the  fire  that  is  actively 


spreading,  then  the  number  of  firelets  corresponds  to  the  perimeter 
or  “front”  of  the  fire.  If  the  birth  rate  and  death  rate  in  our  model 
are  exactly  balanced,  each  firelet  will  on  average  spread  to  one 
other  firelet  before  dying  out.  This  is  consistent  with  spread 
governed  by  a  quasi-linear,  one-dimensional,  perimeter  of 
a  compact,  two-dimensional  fire  shape.  Thus,  a  growth  rate 
/?/<5  =  1  may  be  a  particularly  useful  choice  of  parameters  to 
capture  key  features  of  fire  spread. 

In  this  analysis  we  have  explored  a  relatively  narrow  range  of 
growth  rates  around  /?/< 5=1,  since  this  is  the  regime  that  most 
closely  resembles  wildfire  data.  Significantly  higher  values  result  in 
a  monotonically  larger  fraction  of  fires  burning  all  parcels  while 
much  lower  values  result  in  only  small  events.  Neither  of  these 
extremes  captures  the  relevant  features  seen  in  wildfire  data:  an 
abundance  of  small  events,  and  a  heavy  tail  comprised  of  rare, 
large  events.  However,  if  /?/< 5  is  assumed  to  vary  in  time,  the 
variance  of  this  parameter  can  be  expanded  without  necessarily 
creating  anomalous  distributions.  This  will  be  explored  in  future 
work. 

Coupled  Fire  and  Suppression  Dynamics 

We  include  human  suppression  in  the  dynamics  of  the  fire  by 
adding  a  decision  variable  y,  representing  the  total  resources 
assigned  to  the  fire.  This  modifies  the  death  rate  for  a  system 
with  j  firelets,  so  that  it  depends  on  both  the  natural  extinction 
rate  per  firelet,  <5,  and  the  suppression  rate,  y:  fij  =  y  +  5j.  Note 
that  the  behavior  of  this  model  can  be  fully  characterized  by 
the  initial  size  of  the  fire  (r),  the  scaled  growTh  rate  (/f/<5),  and 
the  scaled  suppression  rate  ( y/S ).  The  scaled  suppression  rate 
captures  the  relationship  between  suppression  due  to  resources 
allocated  and  the  natural  extinction  rate.  As  before,  we  drop  the 
specification  ‘scaled’  and  simply  refer  to  the  suppression  rate.  As 
with  observed  wildfires,  a  larger  fire  requires  more  suppression 
resources. 

Tradeoff  One:  Initial  Size  versus  Suppression  Rate 

Figure  2B  illustrates  P(F)  for  dangerous  conditions  and 
varying  suppression.  When  y  =  0,  the  simulation  produces  a  pure 
power  law  with  exponent  a=l/2  (Figure  2A).  In  comparison, 
increasing  suppression  results  in  monotonically  smaller  fires.  For 
all  of  the  illustrated  values  of  y,  the  statistics  remain  well 
described  by  power  laws,  extending  to  large  event  sizes,  although 


Figure  2.  Fire  size  distributions  P(F)  calculated  for  varying  growth  and  suppression  rates.  Results  are  based  on  104  simulation 
runs  for  a  fire  in  a  region  containing  burnable  substrate  of  7=  1000  parcels.  In  (A),  growth  rate  /1/S  varies  in  the  absence  of  suppression.  Dangerous 
conditions  are  described  by  a  pure  power  law  with  exponent  a=  1/2.  Extreme  conditions  {/!/ 8>  1)  lead  to  excess  weight  in  the  tail,  which  can  be 
interpreted  as  fires  that  have  surpassed  a  threshold  value  and  can  no  longer  be  controlled.  Finally,  mild  conditions  (/J/<5 <  1 )  result  in  a  sharp 
(exponential)  cutoff  Fmax,  which  decreases  as  fl/S  decreases.  In  (B),  the  suppression  rate  y/S  varies  in  dangerous  conditions  (P/S=  1).  Pure  power 
laws  are  obtained  for  a  range  of  y/S,  with  an  increasing  exponent  (steepening  slopes)  as  suppression  is  increased. 
doi:1 0.1 371 /journal. pone.0033285.g002 
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the  exponents  increase  (the  curves  become  steeper)  as  y  increases. 
This  implies  that  the  probability  of  a  large,  system  spanning 
event  remains  non-negligible.  In  contrast,  decreasing  the  natural 
growth  rate  in  Figure  2A  results  in  progressively  smaller 
exponential  cutoffs  in  the  tail  of  P(F),  so  the  probability  of 
extremely  large  fires  becomes  increasingly  remote.  Thus,  varying 
the  suppression  rate  does  not  affect  the  largest  events  as 
significantly  as  varying  the  growth  rate,  reflecting  the  dominance 
of  environmental  factors  over  a  fixed  suppression  allocation  as 
the  number  of  firelets  becomes  large. 

The  distributions  in  Figure  2A  and  2B  are  for  simulated  fires 
initiated  with  j  =  1 .  Here,  the  implicit  assumption  is  that  the  delay 
for  the  arrival  of  suppression  resources  is  negligible.  In  more 
realistic  situations,  the  fire  is  likely  to  have  grown  to  some  larger 
size  by  the  time  suppression  resources  arrive,  and  this  “initial”  size 
(i.e.,  at  the  time  of  arrival)  may  determine  whether  or  not  the  fire  is 
contained.  In  what  follows,  we  vary  the  initial  size  as  a  proxy  for 
fire  growth  during  this  delay,  starting  the  simulations  with  j  =  S 
firelets  and  with  J  —  s  unburned  parcels. 

We  define  the  bum  probability  P(J)  to  be  the  probability  that  the 
fire  consumes  the  entire  burnable  region,  and  investigate  P(J)  as 
a  function  of  the  initial  size  (s),  growth  rate  (yS / <5),  and  suppression 
rate  (y/< 5).  Figure  3  illustrates  our  results  obtained  by  simulating 
100  fires  for  each  set  of  parameters.  Panels  A— C  in  Figure  3 
compare  the  results  for  mild,  dangerous,  and  extreme  conditions. 
In  each  case,  the  larger  the  initial  size,  the  greater  the  suppression 
rate  required  to  contain  it.  Comparing  the  results  for  different 
conditions  shows  that  each  10%  increase  in  the  growth  rate 
requires  roughly  an  order  of  magnitude  decrease  in  initial  size  in 
order  for  the  same  suppression  rate  to  be  effective  in  keeping  the 
burn  probability  low. 

Risk  curves  [46]  generalize  the  concept  of  burn  probability  P(  J) 
to  include  the  full  range  of  sizes  F.  They  are  essential  ingredients 
for  any  interactive,  dynamic  decision  making  tool,  which  might  be 
used  to  estimate  the  suppression  rate  required  to  maintain  the  fire 
at  a  certain  risk  level  based  on  the  current  conditions  and  the  size 
of  the  fire.  Here  P(F)  corresponds  to  the  risk  that  the  fire  will  reach 
size  F  or  greater.  Figure  4  illustrates  the  results  under  mild, 
dangerous,  and  extreme  conditions.  As  expected,  larger  initial  sizes 
s  result  in  increased  risk  of  a  large  fire.  As  the  suppression  rate 
increases  (going  from  solid  to  dashed  to  dotted  lines),  the  risk  of 
large  fires  decreases.  In  mild  conditions  the  initial  size  plays  the 
dominant  role  in  determining  P(J),  and  the  suppression  rate  has 
relatively  little  impact.  In  dangerous  conditions  the  suppression 
rate  has  a  much  more  significant  effect,  especially  for  fires  with 
a  larger  initial  size.  Finally,  in  extreme  conditions  the  risk  curves 
flatten  out,  reflecting  the  likelihood  of  a  size  J  event  for  a  wide 
range  of  initial  sizes  and  suppression  rates. 


Tradeoff  Two:  Time  Delay  versus  Suppression  Level 

In  practice,  the  decision  maker  will  not  know  in  advance  the  size 
of  the  fire  when  suppression  resources  arrive.  This  size  will  depend 
on  the  initial  size  of  the  fire  when  discovered,  prevailing 
environmental  conditions,  and  the  length  of  time  that  the  fire 
grows  before  resources  arrive.  Typically,  a  longer  time  delay 
results  in  a  larger  fire  and  a  need  for  more  resources  to  maintain 
a  low  burn  probability. 

We  evaluate  the  tradeoff  between  the  time  delay  and  resource 
requirements  by  adding  an  explicit  time  dependence  y(t).  For 
simplicity,  we  take  y(t)  to  be  a  threshold  function  in  time.  Letting 
At  denote  the  time  delay,  we  have  y{t)  =  0  for  te[0,At),  and  y(t)  =  y 
for  t>At.  The  function  y{t)  could  be  made  more  complicated,  to 
model  different  resources  arriving  front  different  locations  at 
different  times.  However,  even  this  minimal  form  is  consistent  with 
a  simple  response  policy:  once  a  fire  is  detected,  send  an  initial 
response  force  of  some  standard  size.  In  general,  the  arrival  time 
will  depend  on  the  distance  between  the  fire  and  the  fire  station 
and  the  difficulty  of  traversing  the  terrain. 

Figure  5  shows  the  tradeoffs  between  time  delay  and 
suppression  rate;  we  plot  burn  probability  as  a  function  of  At 
and  y,  under  dangerous  conditions.  Even  for  instantaneous  arrival 
(At  =  0)  P(J)  is  non-negligible  for  small  values  of  the  suppression 
rate  (roughly  y<  1),  suggesting  that  an  initial  response,  no  matter 
how  prompt,  may  be  ineffective  if  it  is  too  small.  For  small  values 
of  At,  contours  of  constant  P(J)  are  roughly  linear  in  semi-log 
space:  log(y)cczl/.  This  indicates  that  the  suppression  strength  y 
must  increase  on  average  exponentially  with  the  delay  A  t  to  keep 
pace  with  the  growing  fire. 

Beyond  a  characteristic  time,  the  slope  of  the  constant  P(J ) 
contours  increases  sharply  with  At,  and  P(J)  becomes  indepen¬ 
dent  of  the  suppression  level.  The  system  enters  this  regime  at 
approximately  a  suppression  rate  of  y  =  100.  This  suppression  level 
roughly  corresponds  to  \[j ,  the  entire  actively  burning  region  (or 
“perimeter”)  of  a  fire  that  burns  all  parcels.  Hence,  at  this  point 
the  suppression  effort  either  reaches  the  fire  in  time  and  is 
extremely  effective,  or  the  fire  has  already  consumed  the  entire 
burnable  region  and  suppression  has  no  effect. 

Tradeoff  Three:  Resource  Division  for  Two  Simultaneous 
Fires 

Dynamic  decisions  involving  resources  are  especially  critical 
when  considering  two  or  more  wildfires  at  different  locations, 
spreading  at  different  rates.  In  many  situations  a  portion  of 
resources  may  be  directed  to  a  new  outbreak,  even  if  they  are 
currently  in  use  elsewhere  on  another,  possibly  larger,  fire.  The 
rationale  is  that  early  response  is  generally  a  priority  because  it 
increases  the  likelihood  of  containment. 
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Figure  3.  Burn  probabilities  P{J)  as  a  function  of  initial  size  sand  suppression  rate  y/< 5.  Conditions  are  mild  in  (A),  dangerous  in  (B), 
and  extreme  in  (C).  As  the  conditions  become  increasingly  severe,  the  transition  from  low  to  high  P(J)  shifts  to  smaller  initial  sizes.  Results  shown 
represent  averages  for  100  fires. 
doi:10.1371/journal.pone.0033285.g003 
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Figure  4.  Risk  curves  P(F)  as  a  function  of  initial  size  sand  suppression  rate.  Conditions  are  mild  in  (A),  dangerous  in  (B),  and  extreme  in 
(C).  In  each  case  the  line  style  indicates  suppression  rate:  y/<5  =  1  (solid),  2  (dashed),  and  3  (dotted).  Burn  probability  P(J)  corresponds  to  P(F  =  J), 
with  J  =  1000.  In  mild  conditions  risk  generally  remains  low,  unless  the  initial  size  is  large.  In  dangerous  conditions  suppression  plays  a  significant  role 
in  reducing  risk  across  the  full  range  of  initial  sizes.  In  extreme  conditions  suppression  is  most  effective  when  the  initial  size  is  small.  Beyond  a  certain 
size,  the  fire  is  increasingly  likely  to  grow  out  of  control. 
doi:1 0.1 371 /journal. pone.0033285.g004 


To  investigate  tradeoffs  involving  the  division  of  resources,  we 
consider  two  fires  (denoted  A  and  B).  We  implicitly  incorporate  the 
notion  of  delays  in  ignition  and  response  in  the  initial  sizes  of  the 
two  fires,  Sa  and  Sg,  taking  Sa  >Sb ,  so  that  Fire  B  may  be  thought 
of  as  the  new  fire.  We  assume  the  conditions  for  Fire  A  are 
dangerous,  and  we  vary  the  conditions  for  Fire  B  from  mild  to 
extreme.  We  place  a  cap  of  y  =  10  on  the  total  resources  available, 
to  be  divided  between  the  two  fires:  y,4+yg  =  y.  Without  loss  of 
generality,  we  set  <5,4=(5g=l.  We  seek  the  optimal  division  of 
suppression  resources  between  these  fires,  where  optimal  is 
determined  by  minimizing  the  combined  size  of  the  two  fires 
(Fa  +  Fb )  (i.e.,  as  a  proxy  for  damage),  averaged  over  1000  runs. 

Figure  6  summarizes  our  results.  Panel  (A)  explicitly  illustrates 
the  optimization  process  when  the  conditions  for  Fire  B  are 
extreme.  The  growth  rate  in  panel  (A)  corresponds  to  the  limiting 
case  (right  most  boundary)  shown  in  the  remaining  panels, 
representing  the  optimal  resource  division  (B),  the  total  burn  area 
of  the  pair  of  fires  (C),  and  the  difference  in  burn  area  between  the 
optimal  solution  and  the  worst  case  scenario  (D). 


Figure  5.  Tradeoffs  between  time  delay  (At)  and  suppression 
rate  (y)  represented  by  contours  of  constant  burn  probability 

P(J).  At  small  time  delays,  y/S  must  exceed  some  minimum  value  to 
keep  P(J)  low.  As  At  increases,  initially  suppression  rate  must  increase 
roughly  exponentially  with  delay  to  maintain  a  constant  P(J).  At  a  larger 
value  of  At,  there  is  a  sharp  transition  where  the  contours  become 
essentially  vertical,  corresponding  to  a  delay,  determined  by  the 
underlying  rates  and  system  size,  beyond  which  the  fire  is  likely  to  have 
either  burned  out  or  grow  out  of  control.  Results  shown  represent 
dangerous  conditions  fS/t 5  =  1  and  averages  over  1000  simulations  and 
all  fires  are  initialized  with  size  s  =  20. 
doi:10.1371/journal.pone.0033285.g005 


For  each  curve  in  Figure  6A,  optimization  identifies  the 
fractional  resource  allocation  yg/y  to  Fire  B  that  minimizes  the 
mean  total  size.  When  Sb  is  small,  the  optimal  solution  is  a  small 
allocation  to  Fire  B,  but  as  .S’/j  increases  the  optimal  allocation  to 
Fire  B  increases  rapidly  and  then  shifts  completely  back  to  Fire  A 
as  sb  becomes  comparable  to  sa-  In  that  case,  because  Fire  B  is 
spreading  more  rapidly,  there  is  a  greater  advantage  in 
concentrating  suppression  resources  on  Fire  A. 

Figure  6B  shows  that  this  general  pattern  persists  while  the 
conditions  for  Fire  B  remain  extreme.  The  contours  shift  upwards 
as  the  severity  of  the  conditions  for  Fire  B  decreases,  implying 
comparable  resource  assignments  become  optimal  at  larger  initial 
sizes  of  Fire  B.  Interestingly,  the  final  transition  from  full  allocation 
to  Fire  B,  back  to  full  allocation  to  Fire  A  is  extremely  abrupt 
(represented  by  the  discrete  color  shift  in  the  upper  right  corner  of 
Figure  6B),  corresponding  to  a  discontinuous  transition  from 
a  solution  that  employs  all  available  resources  to  combat  the 
smaller,  more  rapidly  growing  Fire  B,  to  a  solution  that  gives  up  on 
Fire  B,  and  uses  all  the  resources  for  Fire  A.  When  the  conditions 
for  Fire  B  are  dangerous  or  mild  (Pb/Pa  —  1  )>  the  solutions  are 
similar,  although  in  that  case  it  is  never  advantageous  to  move 
resources  back  to  Fire  A  as  Sb  increases,  because  now  yg  is  the 
more  slowly  growing  fire,  so  when  the  sizes  become  comparable, 
resources  are  most  effective  for  Fire  B. 

Figure  6C  illustrates  the  combined  size  < Fa+Fb )  correspond¬ 
ing  to  the  optimal  resource  divisions  computed  in  6B.  Here  values 
are  strictly  bounded  below  by  the  sum  of  the  initial  sizes  of  the  two 
fires  (recall  .?^  =  100,  and  l<.?g<100),  and  bounded 

above  by  2J  =  2000.  The  transition  from  small  (lower  left)  to  large 
(upper  right)  combined  fire  size  is  a  relatively  smooth  crossover, 
unlike  the  optimal  resource  allocation  contours  in  Figure  6B, 
which  exhibit  sharp  crossovers  and  abrupt  transitions.  In  the  case 
of  two  large  fires,  the  initial  size  and  growth  rate  of  Fire  B  is 
sufficiently  large  that  the  burn  probability  for  Fire  B  approaches 
unity.  In  this  case,  the  size  of  Fire  B  is  typically  near  J,  and  the 
resources  are  concentrated  on  suppressing  Fire  A,  which  on 
average  is  contained  at  a  size  near  3J/5. 

T o  estimate  the  impact  of  optimal  decision  making,  we  compare 
the  combined  size  of  the  optimal  solution  in  Figure  6C  to  the  size 
that  corresponds  to  assigning  resources  in  a  manner  that  maximizes 
the  combined  size  of  the  two  fires  (this  corresponds  to  choosing  the 
value  of  >’g / y  in  Figure  6A  that  maximizes,  rather  than  minimizes, 
mean  total  size  for  each  curve).  This  worst  case  scenario  enables  us  to 
identify  regimes  in  which  optimization  leads  to  the  most  significant 
gains.  Figure  6D  illustrates  the  difference  between  the  maximum 
and  minimum  combined  size.  Interestingly,  the  net  gain  obtained 
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Figure  6.  Resource  optimization  for  two  simultaneous  fires.  Fire  A  and  Fire  B,  subject  to  the  constraint  yA+yB=y=  10.  Fire  A  has 
dangerous  conditions  and  initial  size  ,s^  =  100.  In  panel  (A)  the  optimization  procedure  is  shown  explicitly  for  pB=  1.1  (extreme  case),  and  varying 
values  of  sB.  Optimal  solutions  correspond  to  the  minimum  mean  total  size  (FA  +  F/t>  for  each  curve,  and  begin  with  a  small  allocation  to  Fire  B 
when  sB  is  small,  rapidly  shifting  resources  to  Fire  B  as  sB  increases,  and  then  shifting  resources  back  to  Fire  A  when  sB  becomes  sufficiently  large. 
Results  represent  averages  over  1000  realizations.  In  panel  (B)  these  results  are  extended  to  values  of  fjB  extending  from  mild  to  extreme  conditions. 
The  trends  in  panel  (A)  persist  for  all  cases  where  Fire  B  is  extreme.  When  Fire  B  has  dangerous  or  mild  conditions,  resources  are  not  shifted  back  to 
Fire  A  when  sB  becomes  large.  Panel  (C)  illustrates  the  value  of  (FA+FBy  associated  with  the  optimal  solution  in  panel  (B).  Here  ( FA+FB )  serves  as 
a  proxy  for  total  damage.  Figure  (D)  shows  the  size  difference  between  solutions  corresponding  to  the  maximum  (worst  case)  and  minimum  (optimal) 
value  of  < Fa  +Fb >.  Regions  where  the  difference  is  elevated  correspond  to  regimes  where  optimal  decision  making  has  the  most  significant  impact. 
doi:10.1371/journal.pone.0033285.g006 


through  optimization  exhibits  three  distinct  regions  with  elevated 
values.  One  occurs  along  the  lower  boundary,  where  the  optimal 
solution  allocates  only  minimal  resources  to  Fire  B,  while  the  worst 
case  solution  corresponds  to  assigning  most  of  the  resources  to  Fire 
B.  Another  peak  occurs  in  the  upper  right  corner,  where  the 
optimal  solution  corresponds  to  allocating  all  the  resources  to  Fire 
A,  because  Fire  B  is  sufficiendy  large  and  rapidly  growing  so  that 
the  resources  have  minimal  impact  and  the  burn  probability  is 
essentially  equal  to  one.  Flere  the  worst  case  solution  leaves  all  the 
resources  with  Fire  B,  so  that  in  the  end,  both  fires,  rather  than  just 
one,  tend  to  burn  to  size  J.  The  third,  largest  and  broadest  peak  in 
Figure  6D  occurs  near  the  lower  right  corner,  corresponding  to 
a  small  initial  size  and  rapid  growth  rate  for  Fire  B.  It  is  in  this 
region  where  judicious  allocation  of  resources  has  the  greatest 
impact  on  keeping  both  fires  under  control. 

Discussion 

Recent  natural  disasters  worldwide  highlight  the  urgent  need  for 
real-time,  dynamic,  decision  making  tools  for  response  [47].  In  the 
past,  there  has  been  extensive  work  developing  hazard  maps,  and 
other  static  tools  that  chart  current  conditions  and  estimate  the 
probability  of  regional  hazards  for  policy  makers  and  responders 
[48-50].  In  contrast,  the  next  generation  of  dynamic  tools  will 
have  the  capacity  to  evolve  in  real  time  with  both  the  threat  and 


the  response,  and  use  probabilistic  methods  to  estimate  the 
outcome  based  on  the  combined  disaster  and  decision  making 
trajectories. 

In  principle,  the  information  fed  into  these  tools  could  be  as 
detailed  as  desired  and  involve  a  combination  of  observations  and 
modeling.  In  the  case  of  fire,  it  could  combine  physical  models  of 
fire  spread  and  atmospheric  conditions,  with  high  resolution 
measurements  of  topography,  fuel  conditions,  and  critical  in¬ 
frastructure  [51-53],  For  response,  it  could  include  the  mobility, 
potency,  and  cost  of  different  platforms-aircraft,  helicopters, 
bulldozers,  or  engines.  Each  of  these  inputs  has  seen  significant 
advances  in  recent  years-increases  in  computational  capacity  and 
satellite  observations  have  improved  the  resolution  of  models  and 
observations,  and  new  technologies  for  firefighting  and  detection 
have  improved  the  ability  to  sense  and  respond. 

The  work  presented  in  this  paper  begins  to  address  what  we 
perceive  to  be  a  critical  missing  link  in  these  recent  advances— 
development  of  a  robust  framework  for  combining  and  distilling 
dynamic  information  in  a  manner  that  is  useful  for  policy  and 
decision  making.  Although  we  have  focused  on  wildfire  response  in 
this  paper,  we  believe  that  the  need  for  this  connection  exists  more 
broadly  in  other  situations  involving  coupled  dynamics  of  disaster 
and  response.  The  ultimate  end-product  of  this  work  might  be 
a  tool  that  would  enable  information  from  a  variety  of  sources  to 
be  summarized  in  an  evolving  decision  diagram  with  cost 
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projections  carried  forward  over  a  user-designated  time  interval. 
For  example,  such  a  diagram  might  summarize  model  and 
observation  based  projections,  constrained  by  resource  availability 
and  effectiveness,  while  the  user  considers  deployment  in  different 
locations  at  specific  times.  In  the  context  of  the  work  presented 
here,  our  input  parameters:  the  natural  growth  /?  and  extinction  (5 
rates,  the  initial  size  s,  and  the  suppression  allocation  y  represent 
real-time  inputs  based  on  current  conditions,  models,  and  policy 
decisions,  with  dynamic  outputs,  analogous  to  our  probability 
contours,  highlighting  the  essential  tradeoffs  and  estimating  risks 
over  various  time  scales. 

Of  course,  a  great  deal  must  be  done  before  such  complex 
analysis  can  be  automated.  Current  capabilities  to  observe  and 
model  greatly  exceed  the  capacity  to  distill  and  integrate  results  in 
a  manner  that  is  comprehensible,  efficient,  and  robust  for  policy. 
Any  effort  to  compress  information  is  inherently  delicate  because 
of  the  intrinsic  fragilities  associated  with  cascading  breakdowns 
initiated  by  one  small  failure.  This  motivates  our  approach,  based 
on  deliberately  transparent  models  and  isolated  tradeoffs,  aimed  at 
clarifying  fundamental  issues  that  arise  generally  in  this  class  of 
problems. 

In  these  first  steps  of  developing  a  decision  making  framework 
for  policy,  it  is  less  crucial  to  capture  every  detail  of  wildfire 
dynamics  than  it  is  to  correctly  represent  those  features  that  are 
resolved  in  the  initial,  coarse  representation  that  is  used-in  this 
case  spread  rates  consistent  with  compact  fire  footprints  and  the 
resulting  statistical  distribution  of  sizes.  Such  features  are  expected 
to  persist  as  more  detailed  models  are  incorporated  into  the 
framework.  In  contrast,  a  framework  based  on  a  model  that  does 
not  get  the  simple  things  right  is  likely  to  produce  misleading 
results  from  the  start.  More  detailed  models  of  fire  spread  retain 
compact  shapes  and  statistical  features,  while  incorporating 
additional  features  (e.g.,  topography  and  fuel  type),  and  exposing 
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additional  sensitivities  (e.g.,  winds  and  humidity)  that  are 
important  for  fire  spread. 

This  work  suggests  many  directions  for  future  research- 
integration  with  geospatial  approaches,  increases  in  fire  and 
response  model  resolution,  and  more  accurate,  data-driven 
incorporation  of  physical  parameters.  In  the  case  of  California 
wildfires,  a  dominant  factor  is  wind.  The  largest,  most  destructive 
wildfires  are  typically  associated  with  extreme,  high  wind,  low 
humidity  Santa  Alia  conditions  [54,55].  Fluctuations  in  wind 
speeds  during  these  events  dominate  both  the  dynamics  of  fire 
spread  and  the  effectiveness  of  response.  Even  in  our  abstract 
model,  this  could  be  represented  by  time  dependent  birth  and 
death  rates  for  fire  spread,  to  evaluate  the  optimal  time 
dependence  for  effective  response. 

Current  decision  making  for  wildfire  relies  on  a  combination  of 
protocols  established  in  advance  and  real-time  expert  opinion.  A 
great  deal  of  data  is  available  to  facilitate  these  decisions.  In  some 
cases,  it  may  be  even  too  much,  resulting  in  information  overload. 
The  purpose  of  investigations  such  as  ours  is  not  to  replace  the 
experts,  but  rather  to  provide  them  with  a  sound,  statistical  basis 
for  synthesizing  information,  establishing  protocols,  and  predicting 
system  sensitivities  in  advance. 
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